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Prolate spin-weighted spheroidal harmonics play a key role in black-hole perturbation theory. In 

particular, the highly damped quasinormal resonances of rotating Kerr black holes are closely related 

to the asymptotic eigenvalues of these important functions. We here present a novel and compact 

derivation of the asymptotic eigenvalues of the prolate spin-weighted spheroidal harmonics. Our 

. . . analysis is based on a simple trick which transforms the corresponding spin-weighted spheroidal 

CO \ angular equation into a Schrodinger-like wave equation which is amenable to a standard WKB 

T-H , analysis. Our analytical results for the prolate asymptotic spectrum agree with previous numerical 

^^ ' computations of the eigenvalues which appear in the literature. 
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Dh' I. INTRODUCTION. 

<\ 

r~^ ' The characteristic dynamics of test fields in black-hole spacetimes has been studied extensively since the pioneering 

T— I . work of Regge and Wheeler [l[, see also 0-0| and references therein. An astrophysically realistic model of wave 

dynamics in black-hole spacetimes should involve a non-spherical background geometry with angular momentum. In 

f \ , terms of the Boyer-Lindquist coordinates, the spacetime of a rotating Kerr black hole is described by the line-element 
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where M and a are the mass and angular momentum per unit mass of the black hole, respectively. (We use gravitational 
units in which G = c = 1). Here A = r^ — 2Mr + a? and p = r^ -f a^ cos^ 0. 

In this paper we consider perturbations of the non-spherical Kerr spacetime. The dynamics of a test field ^I* in the 
rotating Kerr spacetime is governed by the well-known Teukolsky master equation [7[. One may decompose the field 



Tj- as 

^ ■ iXl, (+ ^ ft rh\ — J"^4' 
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^im{t, r, e, 0) = e""^5z„(0; au) ,^nm{r)e-''^' , (2) 



where w is the (conserved) frequency of the mode, I is the spheroidal harmonic index, and m is the azimuthal 
harmonic index. The parameter s is called the spin weight of the field, and is given by s = ±2 for gravitational 
perturbations, s = ±1 for electromagnetic perturbations, s = ±i for massless neutrino perturbations, and s — Q for 
scalar perturbations 7]. 

;^ , With the decomposition ([2]), i/" and S obey radial and angular equations, both of confluent Heun type [THIH. 

C^ ■ coupled by a separation constant A{auj). The radial Teukolsky equation is given by [3| 



dr\ dr J 



— a u + 2mauj — A + Aisojr 



^ = , (3) 



where K = {r'^ + a^)aj — am. 

The angular functions S{6; auj) are the spin-weighted spheroidal harmonics which are solutions of the angular 
equation 0-[n| 
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sm^— - -H c^ cos^ 6* - 2cs cos 6* - ^ ^ — + s + AS = 0, (4) 
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where c = auj. The angular functions are required to be regular at the poles = Q and 9 — n. These boundary 
conditions pick out a discrete set of eigenvalues {sAim} labeled by the integers m and I. [In the c — ?> limit these 
angular functions become the familiar spin-weighted spherical harmonics with the corresponding angular eigenvalues 
A = l{l + l)~s{s + l) + 0{auj).] 

The spin- weighted spheroidal harmonics S{9;c) and their corresponding eigenvalues {sAi„i} have attracted much 
attention over the years from both physicists and mathematicians [Ill4l7l |. It is worth emphasizing that in order 



to compute the characteristic resonances of black holes 0, one must first compute the (closely related) angular 
eigenvalues {s^zm} [see Eq. ^]. In particular, in the framework of semi-classical general relativity, it has been 
conjectured that the highly damped resonances may shed light on the quantum properties of black holes [13,[ly|- For 
rotating black holes, these highly damped resonances are characterized by c/ — >■ oo (where c/ = 3c). 

The asymptotic limit c/ — >■ oo of the angular equation ^ was studied in [ll|, [l2l for the s = case. It was found 
that the asymptotic scalar eigenvalues are given by 

oAir,, ^ [2{l ~ \m\) + l]\ci\ + 0{1) . (5) 

As pointed out in [lf|, the analysis of [ll|, [ij] for the spin-0 case is somewhat incomplete. The analysis of [III, [l2| 
requires the knowledge of the number of zeros of the scalar harmonics in the interval [0, tt]. However, as emphasized in 
[16| , the approximated solution found in [111 . Il2l | for the c/ — >■ oo limit is only valid in a region far from the end-points 
9 = 0,7r. Thus, the analysis of [ll|, [l2| cannot rule out the possible existence of additional zeros of the angular 
eigenfunctions near the end-points. The possible omission of such zeros would lead to a wrong asymptotic behavior 
of the prolate (with c = icj) eigenvalues, see [lq| for details. In this respect the analysis of [11, J^] is not complete. 
One of the goals of the present paper is to present a more rigorous proof of the formula ([5]) for the asymptotic scalar 
eigenvalues. 

The asymptotic spectrum of the prolate eigenvalues for the general spin case was first studied in Il3ll. However, as 
pointed out in [la, [20| , the analysis of [l3ll f or the general spin case is fundamentally flawed (see [ig, [20| for details) . 
The numerical results presented in Refs. [16ll20l| for the general spin case provide evidence for the asymptotic behavior 

,Aim ^ l\ci\ + 0{1) , (6) 

where the function 7 depends on the spin-parameter s, the azimuthal harmonic index m, and the spheroidal harmonic 
index I [see Eq. (|2ip below]. The main aim of the present paper is to provide a (simple) analytical proof for the 
prolate formula (jG)) in the general spin case. 

II. A COORDINATE TRANSFORMATION 

It proves useful to introduce the coordinate x defined by [13, UM 

a; = ln(tan(-)) , (7) 

in terms of which the angular equation (j?]) becomes a Schrddinger-like wave equation of the form [22l | 

dx 
where 

U{0) ^{m + s cos Of - sin^ e{c^ cos^ 9 - 2cs cos 9 + s + A) . (9) 

Note that the interval e [0, tt] maps into x € [—00, 00]. The Schrodinger-type angular equation ([S]) is now in a form 
that is amenable to a standard WKB analysis. 

III. THE SPIN-0 (SCALAR) CASE. 

We shall first consider the spin-0 (scalar) case. In this case the effective potential U{9) is in the form of a symmetric 
(invariant under the transformation 9 ^ t: — 9) potential well: in the c/ — > 00 limit it has a local minimum at 

O^min = I with C/(oe,ni„) - -A + 0(1) . (10) 

Regions where U{9) < are characterized by an oscillatory behavior of the wave-function S (the 'classically allowed 
regions'), while regions with U{9) > (the 'classically forbidden regions') are characterized by an exponential behavior 
(evanescent waves). The 'classical turning points' are characterized by C/ = 0. There is a pair {o^ ,o^^} of such 
turning points (with o^^ < o^min < o^^)i which in the c/ — >■ oo limit are located in the near vicinity of o^min: 



2 C/5 = , (8) 
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o^^ = ;±^+o(c7^/^). (11) 



A standard textbook second-order WKB approximation for the bound-state 'energies' of a Schrodinger-like wave 
equation of the form ([5]) yields the weh-known quantization condition [23l - t27| 

1, 



dxy/-U{x) = {N+-)TT ■ iV = 0,l,2,... , (12) 

where x are the turning points [with U{x ) = 0] of the potential well, and A'^ is a non-negative integer. 
Using the relation dx/dO — l/sin^, one can write the WKB condition p^ in the form 



^^/^M = (vv+i)^ . ^^0,1,2,.... (13) 

smfc' z 

The WKB quantization condition (|13|) determines the eigenvalues {A] of the associated spin-weighted spheroidal 
harmonics in the large-c/ limit. The relation so obtained between the eigenvalues and the parameters c, tti, s and N 
is rather complex and involves elliptic integrals. However, given the fact that in the c/ — >■ cx) limit the turning points 
0^^ lie in the vicinity oi 9 = ^ [see Eq. pTj) ]. one can approximate the integral in pT|) by |28j 

Q + 



d9^-cj{9-^r + A=iN+^)n ; A = 0, 1, 2, ... . (14) 

Evaluating the integral in p^ is straightforward, and one finds 

^(A) = (2A + l)|c,| + 0(l) ; A = 0,1,2,... (15) 

for the quantized spectrum. This completes our proof for the prolate asymptotic spectrum in the scalar (s = 0) case 

IV. THE GENERAL SPIN CASE. 

We shall now consider the general spin case. In this case the effective potential U{9) is complex-valued. Its minimum 
is located at s^min = f + ^ + ^{^7^)^ while the two turning points are located at 

sO--l±^ + - + OicJ^/^). (16) 

2 |c/| c/ 

A natural generalization of the WKB analysis to the case of complex- valued potentials is provided in j30l| : the WKB 
quantization rule is given by the standard relation 



^^/^M = (7V+1)^ . ^^0,1,2,..., (17) 

smfc' 2 



which can be approximated near 9 = tt/2 by [26 



y° d9^-c]{9 - |)2 - 2isci{9 - |) + A = (TV + 1 )^ ; A = 0,1,2,.... (18) 

As emphasized in t30|, the integration path between the two complex turning points {s9^ ^s ^^} should be chosen such 
that 



{^^U{9)} = (19) 



along the integration contour [30|. In the c/ — > oo limit this requirement is easily fulfilled by a straight line (parallel 
to the real ^-axis) which connects the two turning points p^ . Substituting 9 = 4>— — (where € K runs from q9~ to 
q9^) into pS)) and neglecting terms of order 0(1), one finds that along the path (IT^ the integral (IT51) can be written 
as 

° d4>^-c]{cp-^Y+A={N+^-)T: ■ A = 0,l,2,... . (20) 



This yields 

^(A^) = (27V + l)|c,|+0(l) ; iV = 0,1,2,... (21) 

for the quantized spectrum. This completes our proof for the prolate asymptotic spectrum in the general spin case 
[29j. It is worth emphasizing that the analytical formula (j2ip for the prolate asymptotic spectrum agrees with the 
numerical results presented in [ifl |20] . 
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